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ABSTRACT 



Aims. In this paper we analyse numerically the propagation and dispersion of acoustic waves in the solar-like sub- 
photosphere with localised non-uniform magnetic field concentrations, mimicking sunspots with various representative 
magnetic field configurations. 

Methods. Numerical simulations of wave propagation through the solar sub-photosphere with a localised magnetic field 
concentration are carried out using SAC, which solves the MHD equations for gravitationally stratified plasma. The 
initial equilibrium density and pressure stratifications are derived from a standard solar model. Acoustic waves are 
generated by a source located at the height approximately corresponding to the visible surface of the Sun. We analyse 
the response of vertical velocity to changes in the interior due to magnetic field at the level corresponding to the visible 
solar surface, by the means of local time-distance helioseismology. 

Results. The results of numerical simulations of acoustic wave propagation and dispersion in the solar sub-photosphere 
with localised magnetic field concentrations of various types are presented. Time-distance diagrams of the vertical 
velocity perturbation at the level corresponding to the visible solar surface show that the magnetic field perturbs and 
scatters acoustic waves and absorbs the acoustic power of the wave packet. For the weakly magnetised case the effect 
of magnetic field is mainly thermodynamic, since the magnetic field changes the temperature stratification. However, 
we observe the signature of slow magnetoacoustic mode, propagating downwards, for the strong magnetic field cases. 

Key words. Sun: helioseismology — Sun: magnetic fields — Sun: oscillations — Sun: photosphere — (Sun:) sunspots 



1. Introduction 

The internal structure of sunspots is still not well known. 
Helioseismological techniques, which analyse the influence 
of internal solar inhomogeneities on sound wave propaga- 
tion and the signatures of such waves at the solar sur- 
face, might be of great help in revealing the invisible, 
■ sub-photospheric solar processes. The ability of forward 
numerical simulations to predict and model a number of 
sola r phenomena i n helio seismology has been sh o wn b y 
e.g. IShelvag et all (I2006L l2007f) : iHanasoge etaH (|2007t ): 
iParchevskv fc Kosovichevl (|2007l ) and others. As magnetic 
fields are, perhaps, the most important property of many 
solar features, a new and rapidly developing field is the 
study of the influence of magnetic fields on acoustic 
wave propagation of solar magnetic field concentrations 
such as sunspots or solar active regions. The appearance 
and importance of slow magnetoacoustic waves has been 
shown in forward MHD simulations in polytropi c mod - 
els by [Crouch fc Callvl (|2QQ3h : iGordovskvv fc JaTnl (|2007t ): 
Moradi et (|2QQ9h . Ray- approximation simulations in a 
more realistic and applicable magnetised model have shown 
a sim il ar behaviour to the a coustic waves (Mo radi fc Callvl , 
2008h. lSn"elvag et aTl (|2007l ) have investigated the influence 
of sub-photospheric flows on acoustic wave propagation us- 
ing foward modelling and demonstrated a discrepancy be- 
tween the actual flow profiles and the flow profiles obtained 
by ray-approximation inversion. The simulations of a wave 
packet, constructed from /-modes, which was carried out 



by I Cameron et al.1 ((2008), have shown a good agreement 
with helioseismological observations of sunspots. Now it is 
timely to perform a full forward magneto-hydrodynamic 
simulation of a wave packet propagating through a non- 
uniform magnetic field region in the solar sub-photosphere 
with a realistic temperature profile. 

In the simulations presented here we consider three dif- 
ferent representative configurations of solar magnetic field. 
Each of the cases has a common feature that the field is spa- 
tially localised, allowing a direct comparison of the travel 
speeds and time differences of the wave propagation be- 
tween the magnetised and non-magnetised solar plasma. 
The representative configurations differ in both magnetic 
field strength and geometry. Their spatial structure affects 
the temperature stratification of the simulated sunspot by 
the magnetic tension. We selected two representative strong 
field configurations with opposite effects on the temper- 
ature in the sunspot: one, where the magnetic-field cur- 
vature is strong and, thus, increases the temperature in 
the magnetic field region; and another, where the magnetic 
field curvature is small, and the temperature is decreased 
in the sunspot. The magnetic configurations we apply are 
in magnetohydrostatic equilibrium with the ambient exter- 
nal plasma. These two-dimensional magnetic fields are self- 
similar and non-potential. 

We present the results according to the following general 
structure. Section 2 briefly describes the numerical tech- 
niques that we have used to carry out the simulations. The 
configurations of the magnetic fields and initial conflgura- 
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tion for the simulations are described in Section 3. The 
source used to generate the acoustic modes in the numeri- 
cal domain is presented in Section 4. Section 5 is devoted 
to (i) the techniques of helioseismological analysis we have 
used, and (ii) the results we obtained. Section 6 concludes. 

2. Simulation model 

The cod e SAC (Sheffie l d Ad vanced Code) has been devel- 
oped by IShelvag et al.l (|2008t ) to carry out numerical stud- 
ies. The code is based on VAC (Versatile Advection Code, 
iToth et "all 1 19981 ). however, it employs artificial diffusivity 
and resistivity in order to stabilise the numerical solutions. 
Also, SAC uses the technique of variable separation to back- 
ground and perturbed componen ts to treat gr avitationalfy 
stratified plasma. According to IShelvag et al.l (|2QQ8l ). if a 
plasma is assumed to be in magnetohydrostatic equilibrium 
given by 
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the system of MHD equations governing arbitrary pertur- 
bations of density, momentum, energy and magnetic field 
is written as follows: 
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where p and p b are the perturbation and background den- 
sity counterparts, v is the total velocity vector, e is the 
total background energy density per unit volume, e is the 
perturbed energy density per unit volume, B5 and B are 
the background and perturbed magnetic field vectors, p tb 
is the total (magnetic + kinetic) background pressure, 7 is 
the adiabatic gas index, g is the external gravitational field 
vector, and p t is the perturbation to the total pressure 
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or, in terms of perturbed energy density per unit volume e, 
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Fig. 1. Sketch of the simulation domain geometry, used in 
the simulations. 



Here p t b denotes the total background pressure 



B 2 

Ptb = Pkb + ~y 



(9) 



which, in terms of background conservative variables, gives 

(10) 
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2)51 
J 2 



(11) 



Eqs. ([2j)- ([TT]) are solved using a fourth-order central dif- 
ference scheme for the spatial derivatives and are advanced 
in time by implementing a fourth order Runge-Kutta nu- 
merical method. The source terms D in the right-hand sides 
of the equations denote the artificial diffusivity and resis- 
tivity terms. The simulation domain is shown in Fig.[TJ The 
2D box is 180 Mm wide and 50 Mm deep, and has a reso- 
lution of 960x1000 grid points; the upper boundary of the 
domain is at the solar surface R = Rq. The boundaries of 
the domain are open. The perturbation source is located 
in the upper-middle (500 km below the upper boundary) 
of the simulation box. The synthetic measurement level is 
located at the solar surface. 



3. Magnetic fields and initial conditions 

For an in itial background model, we adopt the Standard 
Model S ([Christensen-Dalsgaard et all Il996l ). The model 
is then adjusted to have the same temperature stratifica- 
tion as the Standard Model, if the constant adiabatic con- 
stant Ti is taken. According to the Standard Model S, the 
pressure at the solar surface R = Rq is equal to p Q = 
7.61 • 10 4 dyn/cm 2 . This creates an upper limit for the ver- 
tical straight uniform magnetic field in the magnetohydro- 
static equilibrium with non-magnetic external plasma to be 
about B max = \/8ttpq = 1.4 kG. The measured magnetic 
field strength in sunspot umbrae is about 2.5—3.5 kG, which 
suggests the implementation of curved magnetic fields in 
the simulations. In the case of curved magnetic field, mag- 
netic tension balances magnetic pressure, thus increasing 
the upper limit for the equilibrium magnetic field. An ex- 
ample of such balance is a potential magnetic field, where 
the magnetic tension is exactly equal to the magnetic pres- 
sure: such a field does not change background pressure and 
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density equilibrium. However, the potential magnetic fields 
have a disadvantage for numerical modelling. If a potential 
magnetic field is considered, the boundaries of the numeri- 
cal domain should be either fixed or periodic to confine the 
magnetic field and prevent it from strong expansion. 

We use a self-similar non-potential magnetic field config- 
urati on (Schliiter & Tem esvaryl . ll958l : ISchussler fc Rempell . 
120051 : ICameron et all l2QQ8h . which can be obtained from 
the following set of equations: 
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Fig. 2. Magnetic field configuration for the case of weak 

(13) magnetic field (Case A) in the sub-photospheric domain of 
the size of 50 Mm in vertical and 180 Mm in horizontal 
direction. The horizontal (B x ) and vertical (B z ) compo- 
nents of the magnetic field are shown. The field lines are 

(14) overplotted. 



where B$ z describes the decrease of the vertical component 
of magnetic field towards the top of the model, and G is 
the function which defines how the magnetic field opens up 
with height. The magnetic field constructed in this way is 
divergence-free by definition. The equilibrium background 
gas pressure and density are then recalculated using the 
magnetohydrostatic equilibrium condition Eq. (pQ). If the 
magnetic field B5 is prescribed, Eq. (pQ) splits into two in- 
dependent equations for the pressure and density deviations 
from the initial state, caused by the magnetic field. These 
equations are then solved numerically in order to obtain 
the gravitationally stratified plasma model with localised 
magnetic field concentration in magnetohydrostatic equi- 
librium. 

Three characteristic situations, mimicking sunspots 
which differ by the magnetic field strength at the visible so- 
lar surface and curvature of the magnetic field, are chosen 
for helioseismic analysis: weak but strongly-curved mag- 
netic field (B Zy Q = 120 G, Case A), strong but weakly- 
curved magnetic field with B z ^ = 3.5 kG (Case B), and 
strongly-curved strong magnetic field (B Zi @ =3.5 kG, Case 
C). The magnetic field structures for these situations are 
shown in Figs. [20 

The curvature of magnetic field changes the temper- 
ature stratification in the domain. For the first case of 
the weak magnetic field (Case A, Fig. [2]), the temperature 
change is rather small (Fig. [5]). Below the temperature de- 
crease, which is caused by the drop of the kinetic pressure 
in the region where the magnetic field is nearly vertical, a 
temperature increase is noticeable. This increase is caused 
by the pressure rise which is needed to compensate the in- 
crease in magnetic tension. 

The two initial configurations with strong magnetic field 
show these processes in greater detail. For Case B (weakly- 
curved magnetic field, Fig. [2]), the temperature is consider- 
ably decreased beneath the solar surface (Fig. [6]), because 
the magnetic field is nearly vertical at the surface. The ef- 
fect of magnetic tension reveals itself in the strongly-curved 
magnetic field configuration (Case C, Fig. [7j). Here, in this 
latter case, the temperature deviation AT/T is mainly pos- 
itive. 
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Fig. 3. Same as Fig. [2j with a strong magnetic flux (B Zi q = 
3.5 kG) but with weakly-curved magnetic field lines, Case 
B. 

4. Acoustic source 

To generate acoustic waves, we introduce a perturbation 
source described by the expression: 



2tt£ 

A sin — - exp 
^0 



(t-Tif 



exp 



(r - r y 



(15) 

where T =300 s, Ti=600 s, <ji = 100 s, <t 2 =0.1 Mm, r is 
the source location. The source is located in the middle of 
the horizontal layer slightly beneath the solar surface (see 
Fig. [Tj). The amplitude of the source Ao is chosen to be 
sufficiently small, which ensures that convective processes 
will not be initiated in the otherwise convectively unsta- 
ble equilibrium, and that the perturbation is still linear, 
i.e. does not change the background strongly. The source 
generates a temporaly localised wave packet with the dura- 
tion of about 600 s, which has the main frequency of about 
3.33 mHz. 

The acoustic response of the simulation box to the 
source is shown in Fig. [U It is evident from the figure that 
the source generates a whole branch of various solar acous- 
tic modes. The p-modes are visible up to high orders. 
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Fig. 4. Same as Fig. [2j with a strongly-curved, strong 
(B Z ^ Q = 3.5 kG) magnetic field, Case C. 
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Fig. 5. Zoom-in of the temperature difference in the mag- 
netic field region for Case A. The temperature increase is 
caused by magnetic tension. 
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Fig. 6. Same as Fig. [5j with a strong magnetic flux and 
weakly-curved magnetic field lines, Case B. Here, the tem- 
perature decreases at the solar surface. This temperature 
decrease is caused by magnetic pressure. Plasma (3 contours 
with levels (3 = 1, 2, 5 are overplotted. 
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Fig. 7. Same as Fig. [5j with the strongly-curved and strong 
magnetic field, Case C. In this case magnetic tension pre- 
vents evacuation of the magnetic region, and the tempera- 
ture is increased. 

Also, in order to check the validity of the simulations, 
the one-dimensional calculation of the eigenmodes of the 
initial background model was performed. The correspond- 
ing eigenfrequencies are overplotted in the figure (solid 
lines). A small difference between the model and the ob- 
served frequencies is caused by discrepancies between the 
Standard Model S and the model implemented here, which 
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Fig. 8. Power spectrum of the vertical velocity perturba- 
tion generated by the source. The p modes are visible up 
to high orders. Eigenmodes of the background model are 
overplotted by solid lines. 
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Fig. 9. Cross-correlation function deduced from the verti- 
cal velocity perturbation at the solar surface, generated by 
the acoustic source. Three wave bounces are clearly visible. 
Group travel times for the first three bounces deduced from 
ray theory are overplotted. 



is accountable for using the equation of state for an ideal 
gas. However, we leave the development of non-ideal equa- 
tion of state, including the ionisation processes, and correc- 
tion of these discrepancies for a future work, since they will 
not play any major role in the calculations of the influence 
of non-uniform magnetic fields on acoustic wave propaga- 
tion in the solar photosphere. 

Next, Fig.[9]shows the time-distance diagram computed 
at the simulated solar surface via cross-correlating the ver- 
tical velocity component, generated by the acoustic source. 
Three wave bounces are clearly visible on the plot. Some 
weak and artificial reflection from the side boundaries, 
which is caused by not perfectly transparent boundaries 
of the numerical domain, is also noticeable. 
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5. Time-distance analysis 

It has been mentioned above that the acoustic source is 
located in the middle of the horizontal layer close to the 
solar surface. This allows us to study the influence of the 
magnetic field on the acoustic response of the simulated 
solar sub-photosphere by comparing the plasma velocities 
to the left (non-magnetic part of the domain) and to the 
right (where the magnetic field is implemented) from the 
source. Similar technique ha s already been proposed and 
used by IShelvag et al.l (|2QQ7| ) to reveal the discrepancies 
between the real and inverted velocity profiles for the sub- 
photosphere with embedded sub-photospheric horizontal 
flows. 

The vertical velocity differences are computed between 
the points, located at the same distance to the left and to 
the right from the source. The difference images, obtained 
in this way, revail the phase shifts and amplitude changes 
the wave packets experience due to the propagation in the 
magnetised region, compared to the non-magnetised one. 

In each case of a magnetic configuration we compute 
the acoustic power of the vertical component of the velocity 
oscillations over the period of the simulation as function of 
horizontal and depth coordinates, 

a p (x,z) = J v 2 z (x, z,t)dt, 

and then consider the ratio between the corresponding 
points of the quiet Sun and perturbed parts of the model. 
A cut at the surface level corresponds to the acoustic 
power measurements deduced observationally. In addition, 
we have measured the travel time perturbations by cross- 
correlating the velocity signal at the source location with 
the signal at the target location, taking the quiet Sun cross- 
correlation func tion as a reference and usin g both Gabor 
wavelet fitting ([Kos ovichev & Duvall, fl997f) or linearised 
definition outlined by Gizon & Birch (2002). 

In the sections I5.lti5.3l we analyse by these means the 
wave propagation through the three cases of magnetic field 
structures. In the section [5^41 the results are compared. 

5.1. Weak magnetic field (Case A) 

The analysis shows that for the configuration of a weak 
magnetic field (Fig. [2]) the influence the magnetic field exer- 
cises on the wave propagation is mainly caused by temper- 
ature (and hence local sound speed) changes in the magne- 
tised region. The temperature increase below the simulated 
sunspot (see Fig. [5]) causes negative phase shifts of the wave 
packets propagating through the magnetic structure. 

In Fig. [10] small phase shifts are observed in all of the 
bounces, however, the first bounce wave is only affected 
in the magnetic field region (the phase shifts for the first 
bounce at 0-30 Mm and at 70-90 Mm are of the order of 
numerical noise). It was found that higher order bounces are 
all affected from 30 Mm distance onwards from the source. 

In Fig.fTTlwe present the acoustic power ratio measured 
in the vertical velocity as a function of two spatial coor- 
dinates between the waves propagating in the quiet Sun 
and in the perturbed part of the model. We found that the 
power deficit at the location of the flux tube is confined to 
the uppermost layers of the model with the acoustic power 
ratio falling the closer to the surface one measures at. The 
ratio is constant and equal to unity to the left of the flux 
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Fig. 10. Vertical speed difference image for Case A. 
Difference values are computed between points located at 
the same distance but opposite sides of the source. Two 
dashed lines bound the magnetic region with \B\ > 25 G. 
The first bounce (leftmost in the figure) is affected only lo- 
cally by the magnetic field, however, the second and third 
bounces are also affected in the 60-80 Mm distance region. 




Fig. 11. Synthetic acoustic power ratio image for the 
weakly magnetised region of the weakly magnetised model 
(Case A). The ratios are computed between the points lo- 
cated at the same distance and opposite sides from the 
source. The image shows the regions of acoustic power de- 
creased compared with the ambient non-magnetic medium. 
The black lines are the contours of vertical magnetic field 
at 50 and 100 G, respectively. The acoustic rays for e.g. the 
frequency / = 4.5 mHz computed for the quiet Sun model 
are also overplotted. 



tube, while to the right we observe variations in the acous- 
tic power at depths of 3 Mm and lower, extending along 
the straight lines starting below the surface at the flux tube 
boundary. From the overplotted acoustic rays computed for 
the unperturbed model we find that these lines agree well 
with the second and third bounce ray envelopes. 

5.2. Weakly-curved, strong magnetic field (Case B) 

In Case B (see Fig. [3]), the temperature distribution is 
such that the temperature decreases in the sunspot region 
(Fig. [6]), similarly to Case A, since the magnetic tension is 
relatively small. Again, as for Case A, accordingly to the 
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Fig. 12. Same as Fig. [10] but for the strong weakly-curved 
magnetic field (B Zj q = 3.5 kG), case B. Two dashed lines 
bound the magnetic region with \B\> 250 G. 
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Fig. 13. Kinetic energy density ratio for the strong mag- 
netic flux but weakly-curved magnetic field lines, Case 
B. The lower turning points for the first four bounces of 
rays emanating from the source are overplotted in white 
color. The black lines are the contours of the magnetic field 
\B\ = 200 G. 



vertical velocity difference image Fig. [12j the first bounce is 
affected by the magnetic field only in the magnetic field re- 
gion. It is found that the second and higher-order bounces 
carry information about the interaction with magnetic field 
also in the non-magnetic or weakly magnetised sub-surface 
regions. Intuitively, following the temperature structure, a 
delay in arrival time of the wave packet at a distance from 
the source is expected, since the sound speed in the simu- 
lated sunspot is lower than in the non-magnetic surround- 
ing plasma. A more detailed analysis nicely confirms this 
expectation. The first ridge in the image of velocity differ- 
ence is positive, meaning that the wave arrives later at the 
point in the sunspot when compared with the counterpart 
wave that arrives at the same distance in the non-magnetic 
plasma. 

In Case B, the slow magneto- acoustic mode is also ob- 
served in the domain (see Fig. [T4|) . Generally, the slow- 
wave motions follow the magnetic field structure and shape. 
Further, a suppression of oscillations is observed at the sur- 
face in the magnetised region. At the distance of 40 Mm to 
the right from the source, the amplitude ratio of horizontal 
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Mm 

Fig. 14. Snapshot of the horizontal component of the ve- 
locity field in the upper layers of the domain, taken at 
t = 5130 s of the simulation for Case B. The slow mode 
is visible in the magnetised region beneath the solar sur- 
face between x — 20 Mm and x = 70 Mm. The magnetic 
field lines are overplotted. 



velocity oscillations at the surface to the source amplitude 
Aq is about 0.0001, while at the distance of —40 Mm to 
the left from the source it is more than 0.0004 (note that 
the image is overexposed in order to reveal the small am- 
plitude structures). Thus, a significant part of oscillation 
energy transforms into slow magnetoacoustic wave motion, 
which propagates downwards along the magnetic field lines, 
and is taken out from the surface. 

The kinetic energy density ratio plot (Fig. [13]) shows 
the lines of decreased ratio, similar to Case A, with the 
structure immediately below the surface of the tube show- 
ing greater complexity than before, perhaps, due to the 
effect of magnetoacoustic waves. Fig. [13] also shows the 
curves corresponding to lower turning points for the first 
four bounces of the quiet Sun rays. Note that only the sec- 
ond, third and fourth bounce turning points correspond to 
the ray envelopes (see F ig.fTT ]) , representing caustic surfaces 
(jKravtsov fc Orlovl [l993h . 

5.3. Strongly-curved, strong magnetic field (Case C) 

Strong magnetic tension in Case C (Fig. 0]) changes the 
kinetic gas pressure in a way that the temperature (and 
hence sound speed) in the sunspot region increases (see 
Fig. [7j). As expected, the sound speed increase leads to 
a faster wave propagation through the magnetised re- 
gion. Correspondingly, the vertical speed difference image 
(Fig. [T5]) shows a negative sign of the first ridge in the 
first bounce. Similar to the previous cases, the first bounce 
is affected by magnetic field only in the magnetic region 
(20 — 70 Mm distance), and the second and higher order 
bounces are affected everywhere after 20 Mm onwards. 

For this case we measured the oscillatory power in pres- 
sure perturbation, scaled by inverse square root of the ini- 
tial local density which is presented in Fig.[16l The same 
straight line structure, as in Cases A and B, is observed. 

As in the previous case (Fig. [14j), despite the plasma 
(3 = 1 surface does not appear in the domain, the slow 
magnetoacoustic mode of roughly the same amplitude, as 
for the weakly-curved magnetic field, is observed in the hor- 
izontal velocity component (Fig. [T7]) . This surprising fact 
suggests that the conversion of acoustic waves into slow 
magneto- acoustic ones is nearly as efficient as in the case, 
where the plasma (3 = 1 surface is located in the domain. 
However, here the amplitude of horizontal velocity compo- 
nent in the magnetic region is larger, than in the case of 
weakly-curved strong magnetic field. 
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Fig. 15. Same as Fig. [TOj but for the strongly-curved, 
strong magnetic field (Case C). Two dashed lines bound 
the magnetic region with \B\ > 250G. 




Fig. 16. Oscillatory power in pressure perturbation, scaled 
by inverse square root of the initial local density for the 
strongly-curved, strong magnetic field (Case C). 
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Fig. 17. Same as Fig. [14j for the strongly-curved strong 
magnetic field, case C. The slow mode is visible between 
the overplotted magnetic field lines. 



5.4. Overall comparison of the effect of magnetic field 
structures on wave propagation 

In the previous sections we have shown that a magnetic 
field of the same strength acts differently on acoustic waves 
depending on the geometry and curvature of the field. Here 
we summarise the findings by showing the travel time dif- 
ference and wave packet amplitude dependencies for the 
representative three analysed cases. 

We have measured the travel time perturbations for the 
generated waves when travelling through the flux tube us- 
ing both Gabor wavelet and the Gizon-Birch definitions. 
A similar technique was applied bv iThompson fc Zharkovl 



(2008). The so-called Gizon-Birch travel time difference 
plots are shown in Fig. [18] for all three magnetic field cases, 
respectively. Each plot shows that the sign of travel time 
difference, calculated for the first bounce, changes corre- 
spondingly to the sign of the temperature difference caused 
by the magnetic field curvature. The travel time differ- 
ences, obtained for the simulation with the weakly-curved 
strong magnetic field (Case B), are o f the same order of 
those obtained from the observations (iDuvall et all . [l997; 
iHughes et all . 120051 : IZharkov et all . [2007|). Thislact sug- 
gests that the magnetic configuration, used for this simu- 
lation, is close to the magnetic field structures in the real 
sunspots. 

Fig. HH demonstrates the difference in the acoustic 
power absorption and suppression by the three different 
magnetic field configurations at the solar surface. The green 
curve corresponds to Case A. The oscillatory power sup- 
pression in the magnetic field region reaches 30%, however, 
no energy absorption is observed, since the power ratio at 
the distance x = 80 Mm is close to unity. This fact, to- 
gether with the absence of noticeable slow magnetoacous- 
tic mode in the horizontal component of velocity, confirms 
our suggestion that the weak magnetic fields act only as 
temperature and sound speed perturbations for the waves 
propagting through the fields. 

The power ratio for Case C (Fig. [T9], red curve) is rather 
different. A very strong suppression of plasma motions in 
the magnetic field region (up to 90%) is observed. This 
suppression is caused by the strongly increased tempera- 
ture in the simulated sunspot. Also, the energy absorption 
of about 5% is obtained at x — 80 Mm. This confirms 
the partial conversion of the wave packet energy into slow 
magnetoacoustic mode, which propagates downwards and 
removes the energy from the solar surface. 

An even more complicated behavior is demonstrated by 
the acoustic power ratio at the solar surface, calculated for 
Case B (Fig. [T9j black curve). In this case, the absorption 
reaches the value of about 10% at x = 80 Mm. The charac- 
ter of energy suppression in the magnetic field region is also 
different. The first feature here is that the curve is not sym- 
metric with respect to the vertical axis of the magnetic field 
configuration (note that the green curve, corresponding to 
the weak magnetic field (Case A) is completely symmetric 
around the axis, and the red curve is also nearly symmet- 
ric, if the energy absorption is not taken into account). We 
suggest that the wavy structure between x = 35 Mm and 
x = 50 Mm is connected to the conversion of purely acoustic 
wave into slow magnetoacoustic mode. However, the most 
noticeable feature of this curve is the acoustic power in- 
crease of the order of 7-8% at the distance x = 60 Mm. 
This feature is also clearly visible in the two-dimensional 
power image (Fig. [T3|) . and can be compared with acous- 
tic power haloes around sunspots, u sually observed ( see, 
for example, recent ob servations by iHill et all (|200ll ) or 
Na gashima et al.l (|2007f ) using new instrument Hinode, and 
references therein). 

The straight lines emanating at different angles from the 
flux tubes' left boundary to the right edge of the box, visible 
in both oscillatory power plots for vertical velocity compo- 
nent, kinetic energy and pressure perturbation, in our view 
can be explained in terms of the ray theory as the caustic 
surface changes occurring due to the sound-speed inhomo- 
geneity in the x-direction. There are caustics corresponding 
to the envelope to the ray paths for the second and higer- 
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Fig. 18. Travel time difference plots computed for the first 
bounce for weak (case A), weakly curved strong (case B) 
and strongly curved strong (case C) magnetic field cases. 



order bounces: these are approximately co-located with the 
loci of the lower turning points for the second and higer- 
order bounces, illustrated in Figs. [Til E3 UM 

Viewed as the focusing points for the generated waves, 
the caustics can be char acterised by an increase of the 
power in the oscillations ([Kravtsov fc Orlovl Il993l ). Thus, 
the ratio of the power between the two cases can be ex- 
pected to be most pronounced at such surfaces. This power 
increase is clearly observed in the time series (movies) 
of the s imulated wave field, available in online material 
at http://robertus.staff.shef.ac.uk/publications/acoustic/ . 
The magneto-acoustic mode generation and propagation is 
also clearly visible in the movies. 

6. Conclusions 

In this paper we presented the numerical modelling and 
helioseismological analysis of three physically different, lo- 
calised magnetic field concentrations, mimicking sunspots 
in the solar photosphere. The model photosphere is based 
on the solar Standard Model S. The acoustic response of 
this quiet (non-magnetic) Sun model is close to that of the 
real Sun. The implemented magnetic fields for the simula- 
tions are different not only by their strength, but also by 
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20 40 60 80 

x (Mm) 



Fig. 19. Acoustic power absorption by the simulated 
sunspots. The green dash-dotted, black solid and red 
dashed curves correspond to the case A, case B and case 
C magnetic field configurations, respectively. Horizontal 
dashed line represents the power ratio 1. Vertical dash- 
dotted line shows the axis of the magnetic configuration. 



the curvature of the field lines. The curvature of magnetic 
field creates magnetic tension, which consequently changes 
the pressure, density and temperature stratification of the 
equilibrium model. Three representative cases of magnetic 
fields in equilibrium with the external non-magnetic photo- 
spheric plasma are considered: weak magnetic field, strong 
but weakly-curved magnetic field and strongly-curved but 
strong magnetic field models. As a result of different mag- 
netic field geometry, different temperature structures were 
obtained. In the case of weak magnetic field (Case A), the 
temperature deviation from the background is small, how- 
ever, there is a complex structure of temperature decrease 
in the photosphere and temperature increase in the sub- 
photosphere. The two strong magnetic field cases (B and 
C) have the same magnetic field strength at the surface 
(3.5 kG). The case of the weakly-curved field (Case B) is 
characterised by the temperature decrease below the solar 
surface. However, the strongly-curved magnetic field (Case 
C) makes the temperature increased there. 

The spatial structure of the models we used to carry 
out the simulations is such that it allows direct and easy 
comparison of the behaviour of the waves going through the 
non- magnetic plasma with the behaviour of the waves inter- 
acting with the magnetic field region. For that, we imposed 
magnetic field only in one half of the numerical domain, 
leaving the other half unaffected by magnetic field. 

We analysed the three magnetic field cases by the 
means of local time-distance helioseismology. Synthetic 
time-distance, time-distance difference and travel time dif- 
ference dependencies were calculated from the simulations. 
The dependencies show that the main part of effect of 
magnetic field on the acoustic wave is due to the change 
of the temperature structure in the sunspot. However, we 
also show that there is an energy leakage downwards in 
the model due to the wave mode conversion from purely 
acoustic to slow magneto-acoustic wave motion. 

Despite the fact that the results are intrinsically correct 
up to the order of numerical noise amplitude, we acknowl- 
edge the disadvantage of their somewhat limited applicabil- 
ity. The simulations are carried out for a magnetic field and 
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background model, which are essentially two-dimensional. 
Thus, the main applicability limitation of our results con- 
sists in the energy distribution in the acoustic modes which 
is quantitatively (but not qualitatively, if only acoustic 
and magneto-acoustic waves are considered) different from 
the three-dimensional case. Also, the absorption of acous- 
tic waves by a magnetic region in reality may be different 
from the one presented due to the difference in the acoustic 
energy distribution. However, since the sound and Alfven 
speeds, and the other main magnetohydrodynamic param- 
eters are not affected in any way by the dimensionality of 
the problem, the travel times and travel time differences 
are also independent on the dimensionality. 

The two-dimensional magnetic fields used in the 
simulations presented in this paper can be extended 
to the three-dimensional cylindrically symmetric fields. 
However, simulations of acoustic wave propagation through 
three-dimensional magnetic structures requires significantly 
larger computing resources, so we leave such analysis for the 
nearest future. 
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